Pentagonal nanowires: a first-principles study of atomic and electronic structure 
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We performed an extensive first-principles study of nanowires in various pentagonal structures by 
using pseudopotential plane wave method within the density functional theory. Our results show 
that nanowires of different types of elements, such as alkali, simple, transition and noble metals 
and inert gas atoms, have a stable structure made from staggered pentagons with a linear chain 
perpendicular to the planes of the pentagons and passing through their centers. This structure 
exhibits bond angles close to those in the icosahedral structure. However, silicon is found to be 
energetically more favorable in the eclipsed pentagonal structure. These quasi one dimensional 
pentagonal nanowires have higher cohesive energies than many other one dimensional structures and 
hence may be realized experimentally. The effect of magnetic state are examined by spin-polarized 
calculations. The origin of the stability are discussed by examining optimized structural parameters, 
charge density and electronic band structure, and by using analysis based on the empirical Lennard- 
Jones type interaction. Electronic band structure of pentagonal wires of different elements are 
discussed and their effects on quantum ballistic conductance are mentioned. It is found that the 
pentagonal wire of silicon exhibits metallic band structure. 



PACS numbers: 68.65.-k, 73.63.-b,61.46.+w, 73.90.+f 

I. INTRODUCTION 

Very thin metal wires produced by the tip retracting 
from nanoindentation in scanning tunneling microscopy 
(STM) or by mechanically controllable break junction 
(MCB J) have been subject ^.arjiumber of experimental 
and theoretical studies .ElflBoffiffl In particular, the step- 
wise behavior of the conductance measured in the course 
of wire stretching at room temperature has attracted 
the interests in various fundamental features of quan- 
tum theory, such as the quantization of ballistic electron 
transport in very thin and one dimensional conductors as 
well as Anderson's localization in very long metal wiresa. 
Recorded values of conductance just before the breaking 
of the wire were in the range of the quantum of conduc- 
tance, G = 2e 2 /h. This implies that the smallest cross 
sections of the wire are of atomic dimensions. In fact, the 
conductance of suspended single atom gold wires, which 
have been jjpaduced recently, is measured, to be very 
close to G clll3. As pointed out earlier JidOlij force and 
conductance variations measured concomittantly during 
stretching have indicated a close connection between the 
atomic structure and the step-wise behavior of conduc- 
tance. It is now understood that a complex interplay 
between the quantization of electronic states with level 
spacing larger than room temperature, and pthc stable 
structure having well defined number of atomsli 4 ! and also 
dynamic self-consistent potential in presence of a current 
flow results in the observed step- wise behavior of conduc- 
tance G as a function of stretching. 

Apart from being a potential nanodevice with multiple 



nanodevices, nanowires are important because of their 
exotic and stable atomic structures occurring in differ- 
ent sizes of cross sections. Therefore, a lot of effort is 
being devoted to the production of nanowires that are 
conducting and stable. Metals crystallize in bulk three 
dimensional structures because that is the most stable 
form. If we wish to create one dimensional (ID) systems, 
there is clearly some struggle against nature. Further- 
more, ID periodic metals can suffer Peierls distortion and 
become non-metallic. For finite nanowires, this tendency 
may be suppressed. However, for longer nanowires we 
could end up with ID systems that are either unstable or 
insulating, both undesirable. Therefore, we should con- 
sider structural arrangements, ID or perhaps quasi-lD , 
which have cohesive energies as close to the bulk as pos- 
sible. Our search has led us to pentagonal nanowires, a 
quasi- ID system, where a pair of pentagons sandwich a 
single atom in a local icosahedral structure. The struc- 
ture looks like a pedestal lamp with a pentagonal base, 
is infinitely repeated along the direction perpendicular 
to the planes of the pentagons, is stable and does not 
suffer from Peierls distortion. We performed extensive 
calculations and found that these pentagonal quasi- ID 
nanowires have higher cohesive energies than many other 
ID structural arrangements. 

The pentagonal structure is incompatible with trans- 
lational symmetry, and hence it is not normally seen in 
2D and 3D crystal structures. Strong evidence for five- 
fold symmetric structures appeared in the first-principles 
molecular dynamics simulations where the observatiopa 
of a 13 atom stable icosahedron of Na was reported.113 
The structure can be viewed as a "tiny pentagonal 
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Na atom present on either side of the pentagons. The 
two pentagons share an apex Na atom and hence the 
Nai3 cluster. Several composite structures with pentag- 
onal motifs have also been observed in simulated, anneal- 
ing study of ultra thin Al and Pb nanowiresiiS Subse- 
quently, suspended monatomic chains, strands and he-, 
lical structures have been realized experimentallylJiiJ'Lil 
As the following discussion shows, part of the reason for 
the stability of the periodic pentagonal structures is that 
among several small planar clusters, made of particles in- 
teracting through a two-body Lennard-Jones potential, 
the pentagonal structure has the highest binding energy 
per particle. 

Different regular atomic structures occurring in dif- 
ferent sizes are now a focus of interest of experimental 
and theoretical studies seeking m ore fundamental , up- , 
derstanding of all these structures MBMSN 
Whether the pentagjmainStructures predicted earlier by 
empirical methods! 1 ^r^r 8 !! 29 ! are common to other ele- 
ments and can be understood from more fundamental 
principles have become an important issue. In this pa- 
per we address this question by using the first-principles 
plane wave calculations within the density functional the- 
ory. We carry out state-of-the-art total energy calcula- 
tions for Na, Al, Cu, Au, Fe, Ni, Pb, Si and Xe in two 
different pentagonal structures (four structures for Au), 
and find that the staggered pentagonal structure is a sta- 
ble structure for these elements except for Si. Further- 
more, we compare the energetics with other linear struc- 
tures and perform an extensive analysis of the electronic 
structure and charge density to reveal the origin of stabil- 
ity and electronic properties of the pentagonal structure. 
Finally, we mention the effect of the pentagonal structure 
on the ballistic conductance. 



II. DESCRIPTION OF THE METHOD AND 
ATOMIC STRUCTURE 

First-principles plane wave calculations are performed 
within the supercell geometry using a tetragonal unit cell. 
The axis of the wire is taken along z-axis, and the lattice 
parameter of the wire coincides with the lattice parame- 
ter c of the tetragonal supercell. The lattice parameters 
of the tetragonal cell in the x-y plane are set as a = b = 15 
A so that the interaction between a wire and its peri- 
odic images are negligible. Bloch states are expressed by 
the linear combinations of plane waves with the cutoff 
energy |k + G| 2 always larger than the optimum cutoff 
energy suggested for the ionic ultrasoft pseudopotcntialE2l 
of the clement under study. The Brillouin zone (BZ) i ( 
tegration is performed within Monkhorst-Pack schern 
using (1 x 1 x 20) k points. Results are obtained by gener- 
alized gradient approximationEj (GGA). Preconditioned 
conjugate gradient (CG) method is used for wave func- 
tion optimization. To find the correct ground state we 
also performed spin-polarized calculations for nanowires 



both VASPE3 and CASTEPEa codes independently. 

We considered the following pentagonal structures: (i) 
Atoms form parallel pentagons which are perpendicu- 
lar to the (z-) axis of the wire with separation w, but 
successive pentagons are rotated by ir/5. In addition, 
a monatomic chain along the z-axis passes through the 
center of pentagons, where each chain atom is located at 
a point equidistant from the planes of pentagons. The 
lattice parameter c in the direction of the chain is twice 
the spacing between pentagons, i.e. c — 2w. This struc- 
ture is specified as the staggered pentagon, S, and has 12 
atoms in the unit cell, (ii) Same as (i), but successive 
pentagons are not rotated so that they have the same 
orientation relative to the z-axis, and hence c = w. This 
structure is the eclipsed pentagon, £, with 6 atoms per 
unit cell, (iii) Staggered pentagonal structure without 
the monatomic chain passing through the centers of the 
parallel pentagons, This is called the 1Z structure and 
has c = 2w. (iv) We also found a modified version of the 
iS structure in gold nanowires, which is specified as the 
deformed staggered pentagon, T>S. Here, adjacent pen- 
tagons are staggered, but one of the two atoms of the 
central chain in a unit cell is slightly displaced, while 
the other one is missing. Accordingly, 1Z and T>S struc- 
tures have 10 and 11 atoms in a unit cell, respectively. 
The pentagonal structures and their relevant structural 
parameters are schematically described in Fig. [j]. 



III. RESULTS AND DISCUSSIONS 
A. Optimized structures and energetics 

The energetics and atomic structure of Na, Al, Pb, 
Cu, Au, Ni, Fe, Si and Xe wires in S and £ have been 
investigated. Our results for the structural parameters 
and binding energies are listed in Table |. The binding 
energies and the relevant inter-atomic distances corre- 
sponding to equilibrium bulk crystal structures are also 
given for the sake of comparison. The binding energy 
per atom for a given structure is calculated as the dif- 
ference of the energy E a of an individual atom and the 
total energy of wire Et having n atoms in the supercell 
divided by n, i.e. Eb — E a — Ex/n. In spite of the fact 
that Na, Al, Cu, Pb, Au, Fe, Ni, Si, and Xe atoms have 
different electronic configurations and form bulk crystals 
with dramatically different properties, they all form sta- 
ble wires in the pentagonal structure. It is very inter- 
esting that Al with 3s 2 3p valence states and Na with 
3s valence state form similar pentagonal structures. Sil- 
icon, a group IV element which is normally crystallized 
in (tetrahedrally bonded) diamond structure is predicted 
to form pentagonal wires similar to what Xe, having a 
closed shell structure, does. The pentagonal structure is 
a stable structure corresponding to a local minimum on 
the Born-Oppcnheimer surface. 

Calculated binding energies show that among the pen- 
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TABLE I: Comparison of calculated structural parameters and binding energy, Eb, for different pentagonal structures of 
different elements. The nearest neighbor distance, d , and binding energy, Eo, are calculated for the optimized bulk crystals. S, 
E, R, DS, and SM are staggered, eclipsed, staggered without central chain, deformed staggered, magnetic staggered, respectively. 
Bond lengths and energies are in A and eV respectively. 
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getically more favorable for Na, Al, Cu, Pb, Au, Fe, Ni 
and Xe. However, the differences in binding energies, 
AEb = Eb,s ^ Ebx are generally small and are in the 
range of ~ 10 meV. AEb < for Si, which favors the 
eclipsed pentagonal structure. 

We investigated the rotation between S and £ struc- 
tures of gold nanowires by breaking the rotation from S 
to £ in seven steps; at each step the relative angle ((p) 
between two pentagons in a unit cell increased by an an- 
gle Aip — 6°. For < (p < 36° the lattice parameter c is 
twice the spacing between pentagons w, i.e. c = 2w. The 
variation of the energy as a function of the rotation angle 
ip is illustrated in Fig. |^. The maximum of the energy 
curve corresponds to Qs^s ■ For robust rotation of pen- 
tagons, the wire is apparently under high compression 
due to the core-core repulsion between two adjacent pen- 
tagons, and consequently Qs^s acquires a high value ~ 
3 eV. However, upon optimizing the structure after each 
rotation step, the structure is modified, in particular, the 
spacing between adjacent pentagons increases with ip. As 
a result, the energies calculated arc lowered dramatically, 
the curve for the variation of the energy as a function of 



cell (or 32 meV/atom). We note that for 30° < <p < 42° 
around the £ structure the energy curve is practically 
flat. This implies that £ structure may acquire a helicity 
along the axis of the wire, if each pentagon rotates by 
a small angle ip. In fact, the helicity in nanowires was 
seen in classical MD calculations!.} and has recently been 
observed experimentallytll. Nevertheless, it is clear that 
the £ structure is not only energetically less favorable, 
but is also unstable. 

From Table | it is easily seen that this is the general 
trend in all systems that show metallic bonding (Na, Al, 
Cu, Pb, Au, Ni). Si also shows a similar behavior, but 
the relaxation of the z axis lattice constant in rotating 
the nanowire from the S to the £ structure is much less 
compared to Na, Cu or Au. This is because, in a system 
with metallic bonding, the electrons are largely delocal- 
ized and they screen the ion cores less effectively. Hence, 
when the cores of the atoms in two successive pentagons 
come closer on the structure being rotated from S to £, 
there is strong repulsion which tend to increase w and 
make the latter structure less favorable. For a system 
like Si, showing directional bonds, there is much better 
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the pentagon forms four bonds-two bonds of 2.57 A with 
two central chain atoms and two bonds of 2.60 A with 
two other atoms in the pentagon. The third neighbor of a 
Si atom is 2.96 A apart, and there is no bond formation 
with it. On the other hand, for a Si atom in the pen- 
tagon for the £ structure, there are 6 bonds-two bonds 
of 2.56 A with two central chain atoms, two bonds of 
2.55 A with two atoms in the pentagon and two weaker 
bonds of 2.73 A with two Si atoms in the pentagons 
above and below it. 

In course of structure optimization of gold wire, we 
found two other structures-7?. and T>S structures. (See 
Fig. [l].) These structures have smaller numbers of atoms 
in their unit cells, and are found to be energetically more 
favorable than the S structure in case of gold nanowire. 
Interestingly, 1Z is found to be energetically less favor- 
able than S for Al nanowires. The stability of the S 
structure is examined for transition metals Ni and Fe by 
performing both spin unpolarized and spin polarized cal- 
culations. We found that Ni and Fe nanowires in stag- 
gered pentagon structures are stable for both nonmag- 
netic (spin unpolarized) and magnetic (spin polarized) 
states. Although, for Ni, the energy gain in the magnetic 
structure compared to the nonmagnetic one is small, the 
spin polarized state of Fe staggered pentagonal wire in- 
creases the binding energy by 0.73 eV/atom. Similarly, 
the magnetic moment per atom is also much larger in 
case of Fe (~ 3/Ltp) than in Ni (~ 0.77/zp). 

The interatomic distance from the chain atom to the 
corners of a pentagon, dc-p, is slightly smaller than the 
nearest neighbor distance within the pentagon, e?p_p, 
but both distances are close to the nearest neighbor dis- 
tance do of the bulk crystal. Assuming that dc-p ~ 
dp-p, the coordination number for a chain atom is equal 
to 10, and that for an atom of the pentagon is 4. Hence, 
roughly speaking, the average coordination number is 7. 
This is smaller than the bulk fee and bee coordination 
numbers 12 and 8, respectively. We note, however, the 
nearest neighbor distance between adjacent pentagons, 
c?p 1 _p 2 ~ 14-20% larger than dc-p and <ip_p. The bond 
angles, ai - as, are close to the bond angles of icosahe- 
dral structure, i.e. 63.4° and 116.6°. Therefore, the local 
atomic configuration piftc pentagonal wires mimics the 
icosahedral structure.t§E3 

Calculated binding energies Ep of pentagonal wires are 
lower than the calculated bulk binding energies Eq in the 
last column of Table |[ This can be explained by higher 
coordination number in bulk crystals. On the other hand, 
the binding energies of pentagonal structure are higher 
than the binding en erg ies of various monatomic chain 
structures. In Table O the binding energies of S struc- 
ture are compared with the binding energies of relevant 
monatomic chain structures calculated^ earlier for Al 
and Au. For example, three monatomic chain structures 
of Al are linear (L), zigzag (W) and equilateral trian- 
gular (T) chains have binding energies 1.87, 1.92 and 
2.50 eV, respectively. The coordination number of the 



TABLE II: Comparison of nearest neighbor distance and 
binding energies of monatomic linear (L), zigzag (W), tri- 
angular (T) chain structures and bulk crystal (B) with the 
binding energy of the staggered pentagon S structure calcu- 
lated for Al and Au. The binding energies of L, W and T 
structures are taken from ReffE4l. 
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ID monatomic chain structures and bulk crystal, so its 
binding energy, Ep=3.2 eV. Gold also follows the same 
ordering. We note the general trend that the binding 
energy increase with increasing coordination number in 
different structures. This trend is clearly observed in Ta- 
ble H by going from the monatomic linear chain to the 
bulk. However, in the Xe pentagonal wires, the bind- 
ing energy becomes larger than the binding energy of 
the bulk fee structure (excluding the contribution of the 
Van der Waals interaction). On the other hand, all the 
nearest neighbor distances (dc-C, dc-p, dp-p, etc.) are 
much smaller than the bulk nearest neighbor distance, 
da- 



B. Energy band structure and total charge density 

The electronic energy band structure of Na, Al, Au and 
Si are presented for both <S and £ structures in Fig. [|. 
Overall forms of the energy bands are the same in both 
structures for Na, Al and Au, except for some shifts and 
splittings of degenerate bands. All these structures have 
the same number of atoms in their cells, but Fermi level is 
crossed by different number of energy bands for different 
elements. Number of bands which cross the Fermi level 
is crucial for the quantum ballistic conductance and the 
stability of nanowire. Under ideal conditions, the con- 
ductance is determined by the number of bands crossing 
the Fermi level, being Go per band. We found that the 
stable staggered pentagon structures of Na, Al and Au 
nanowires have 6, 10 and 6 bands crossing the Fermi 
level, respectively. On the other hand, the Fermi level 
of the unstable eclipsed pentagon structure of Na, Al 
and Au nanowires is crossed by 6, 8 and 5 bands, re- 
spectively. While the staggered pentagonal structure of 
Si nanowire have 6 bands crossing the Fermi level, this 
number is raised to 10 for the eclipsed pentagon, which 
is the stable structure. One notices that a degenerate 
band has dropped below the Fermi energy in going from 
the S to the £ structure for Si. This is possibly because 
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pentagonal planes, which is absent in the S structure, 
and is responsible for the £ structure being more sta- 
ble. While stretching the nanowire, the number of atoms 
in the neck region, where the wire is thinned, and thjd* 
structure exhibit sequential and step-wise changes. EjEZI 
It has been argued that these changes are closely related 
with a band moving jttp from the Fermi level and be- 
coming unoccupied. An important feature of Fig. || 
is that the Si nanowire is metallic in both structures. 
Several bands crossing the Fermi level gives rise to high 
density of states at Ep. We found that all pentagonal 
nanowires except Xe studied in this paper are metallic. 
The Xe nanowire in S structure is a semiconductor with 
a wide band gap. 

The character of the bonding in pentagonal structures 
are revealed by the analysis of electronic charge density. 
In Fig. [| we show the charge density contour plots of Na, 
Al, Au, and Si in different planes. The lateral plane pass- 
ing through the plane of the pentagon show the character 
of the bonding between the atoms in a pentagon. The 
vertical plane includes the central chain as well as one 
atom of the pentagon. Because of inversion symmetry 
of the S structure, one atom of the adjacent pentagon is 
also included in the same vertical plane. The charge dis- 
tribution of Na is uniform in vertical and lateral planes. 
Because of low charge density significant structure cannot 
be resolved. The charge density contour plots in Fig. |] 
display some differences in the charge distributions be- 
tween different atoms. The charge distribution of the 
central chain foxJU is reminiscent of that of monatomic 
chain structured and has a directional character. The 
directional behavior is, however, less pronounced in the 
bonds forming between two atoms in the same pentagon 
and also between central chain atom and pentagon atom. 
The charge distribution of gold wire reflects the charge 
distribution of bulk metal. The directionality with high 
density along the line connecting two nearest neighbor 
atoms is absent. Non uniform charge distribution with a 
directionality between nearest neighbor atoms is clearly 
seen in the contour plots of Si both in S and £ structures. 
The directionality of charge distribution originates from 
the valence of the element which makes the pentagonal 
wire. Here, Al and Si with valence states consisting of 
3s, 3p x ,y,z orbitals form directional bonds. In contrast, 
Na and Au with s valence orbitals exhibit bulk like, uni- 
form charge distribution which is characteristic of metals. 
Despite these differences, all the pentagonal nanowires in 
Fig. H are metals with finite density of states at the Fermi 
level. 



C. Discussions 

As corroborated by the present first principles calcu- 
lations, the pentagonal structure is one of the energet- 
ically favorable structures of wires having translational 



periodicity along its axis. Of course, it is only a lo- 
cal minimum on the Born-Oppenheimer surface and oc- 
curs for a given number of atoms in the cross section 
of the wire (or in the ID unit cell). Since the pentago- 
nal structure occurs for a number of different elements 
as demonstrated in this study, the stability must stem 
from the pentagonal geometry. We examined the rel- 
ative stability of the pentagonal geometry by perform- 
ing a simple analysis based on the two-body Lennard- 
Jones potential, V(r^j) = Ae[(a/r h j) 12 ~ (c/r^) 6 ]. We 
considered three structures enco unte red in the structure 
optimization of very thin wiresjl3t3 namely equilateral 
triangle, pentagon and hexagon. The total potential, 
Vt = J2i^jVi,j( r i,j) f° r these structures are obtained 
by adding all the two-body potential energies, Vij. By 
optimizing the total potential energy Vt relative to the 
nearest neighbor distance d\, dVp/ddi = 0, the opti- 
mized interatomic distances and the binding energies for 
these three structures are calculated in terms of the pa- 
rameters a and e. We found binding energies, Eb of — e, 
-1.1111 e, and -1.09 e, for equilateral triangle, pentagon 
and hexagons, respectively. The corresponding optimized 
nearest neighbor distances are 2 1 / 6 cr, 1.89(7 and 1.91cr, re- 
spectively. It is very interesting to note that among these 
structures, the pentagon has the highest binding energy. 



IV. CONCLUSIONS 

We carried out an extensive analysis of the energetics 
and structure of thin wires made from different elements, 
such as alkali, simple, noble and transition metals, and 
also Si, Xe. First-principles calculations with fully opti- 
mized structures yield that the ID structure formed by 
parallel but staggered pentagons and an atomic chain 
passing through the center of pentagons is generally sta- 
ble and energetically favorable relative to other pentag- 
onal structures. However, there are exceptions. For ex- 
ample, while the eclipsed pentagonal structure is favored 
by Si nanowires, in gold wires, different versions of pen- 
tagonal structures are found to be energetically more fa- 
vorable. The binding energies are intermediate between 
ID chain structure and bulk crystal. All nanowires of 
different elements studied in this paper, except Xe, are 
metallic in the pentagonal structure. Strong cohesion 
and metallicity of quasi-lD pentagonal nanowires sug- 
gest that they can be useful in practical applications and 
deserve further experimental studies. 
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Figure Captions 




FIG. 1: Schematic description of various pentagonal struc- 
tures with the structural parameters: Lattice parameter along 
the wire, c, and spacing between adjacent pentagons, w. (a) 
5 staggered pentagon structure with c = 2w. Numerals spec- 
ify atoms. Atoms 1 and 2 form the chain passing through the 
centers of the staggered pentagons. Relevant inter-atomic dis- 
tances dc-c = w, dc-p, dp~p and dp 1 -p 2 are between the 
atoms (1-2), (1-3), (3-4) and (4-6), respectively. Bond angles 
ai, Q2, Q3, en and as occur between (3-1-4), (3-1-5), (3-1-6), 
(3-1-7) and (3-6-4). (b) In the eclipsed pentagon structure £ 
all the pentagons are aligned. Hence dc-c = dp 1 -p 2 — w and 
c = w (c) Staggered pentagonal structure, 1Z, is similar to 5 
structure in (a), except that central atomic chain is missing, 
c = 2w (d) The deformed staggered pentagon structure T>S 
with c — 2w. 
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FIG. 2: Variation of the total energy with the rotation angle 
ip from S structure (ip — 0) to £ structure (tp — 36°). The 
maximum energy is the energy difference, Qs-t£, between the 
S and the £ structures. Energies indicated by dots correspond 
to the robust rotation of pentagons without relaxation of the 
structure. Energies indicated by diamonds are calculated by 
fully optimizing the structure including the lattice parameters 
at each step. The inset shows variation of c. 
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FIG. 3: Energy band structures of S and £ structures, re- 
spectively, (a) and (b) for Na, (c) and (d) for Al, (e) and (f) 
for Au, and (g) and (h) for Si. Fermi level shown by dashed 
lines marks the zero of energy. 



a) o 



b) 




FIG. 4: Left panels are charge density contour-plots in lateral 
planes which coincide with the plane of pentagons. Right 
panels are in vertical planes which pass through the central 
chain and one atom of each pentagon, (a) and (b): Na in S 
structure, (c) and (d): Al in <S structure, (e) and (f): Au in 
S structure, (g) and (h) Si in S structure, (i) and (j): Si in 
£ structure. 



